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The magnetic and transport properties of the metal phthalocyanine (MPc) and FieMPc (M = 
Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn and Ag) families of molecules in contact with S-Au wires 
are investigated by density functional theory within the local density approximation, including local 
electronic correlations on the central metal atom. The magnetic moments are found to be consider¬ 
ably modified under fluorination. In addition, they do not depend exclusively on the configuration of 
the outer electronic shell of the central metal atom (as in isolated MPc and FieMPc) but also on the 
interaction with the leads. Good agreement between the calculated conductance and experimental 
results is obtained. For M = Ag, a high spin filter efficiency and conductance is observed, giving 
rise to a potentially high sensitivity for chemical sensor applications. 
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I. INTRODUCTION 

Metal pht halo cyanines (MPcs) constitute a family of 
medium sized molecular semiconductors, which are of 
considerable interest for numerous applications such as 
chemical sensors [I|, fuel cells [2|, solar cells Q, and op¬ 
toelectronic devices Q. A large number of MPcs can be 
synthesized Q and geometrically arranged in large areas 
at low cost 0, which is important from the technolog¬ 
ical point of view. In particular, the photoconductivity 
of MPcs has been studied intensively with the purpose 
of increasing the electrical conductivity of devices based 
on these molecules Q. The improvement of the organic 
semiconductor devices relies on the quality of the metal- 
organic interfaces [H, in particular, on the efficiency of 
chame injection and on the mobility of the charge carri¬ 
ers [ 3 , [i3- For example, there is a charge transfer shift of 
the electronic levels at the CuPc/Au interface in the early 
stages of CuPc deposition on An llE . and new occupied 
molecular states are created [l2, [i3 • 

The replacement of the M atom as well as the substi¬ 
tution of H by F (fluorinated MPc) alters the gap be¬ 
tween the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital M , the lat¬ 
ter resulting in n-type conduction, in contrast to the p- 
type behavior of the standard MPcs [il,®. The charge 
transport through the molecule can be measured directly 
1^ . fig or indirectly by photo-induced electron transfer 
, time resolved microwave conductivity [2^ , and scan¬ 
ning tunneling microscopy 0. CuPc has been contacted 
with An atomic chains on NiAl substrate by Nazin et al. 
[22[, who observed a shift and splitting of the molecular 
orbitals as well as modifications of the electrode orbitals 
by scanning tunneling microscopy. Field-effect transis¬ 
tors and metal-insulator-semiconductor diodes have been 


used to study the transport through CuPc for different 
leads such as Ca, An, and F4TCNQ/AU, demonstrat¬ 
ing both electron or hole transport with a strong de¬ 
pendence on the geometry of the molecule-metal con¬ 
tact [2^. CuPc sandwiched between two semi-infinite 
An eiectrodes has been investigated theoreticaiiy in Refs. 
[23, [ill. The transmission coefficient, T{E), shows two 
peaks near the Fermi energy (Ep) which have been dis¬ 
sected in terms of moiecuiar orbitais. The eiectronic 
states of CuPc hardiy change when ieads are attached. 

On the other hand, the Landauer approach and 
Green’s function formaiism have been used to address 
the quantum transport in MPc structures. MPcs with M 
= Mn, Fe, Co, Ni, Cu and Zn sandwiched between semi¬ 
infinite armchair singie-waii carbon nanotubes (SWNT) 
have been considered in Ref. [ 26 | using the SMEAGOL 
package [i^, Within the generaiized gradient ap¬ 

proximation (GCA) it has been conciuded that EePc and 
MnPc can be used as spin fiiters, and that the spin fiiter 
efficiency increases when N-doped graphene nanoribbons 
are used as ieads for EePc [^. To overcome the weak 
interaction between An contacts and MPcs, S atoms have 
been added, ieading to a distinct moiecuiar bonding. The 
transport properties of such AuS-MPc-SAu (M = Cu, 
Mn) systems have been studied by the WanT code, with 
the resuit that eiectronic correiations are iikeiy to be ir- 
reievant [s^. Additionaiiy, the authors have conciuded 
that CuPc is a moiecuiar conductor, and MnPc a spin 
fiiter. 

In this articie we extend previous works on AuS-MPc- 
SAu and AuS-EieMPc-SAu junctions by inciuding spin 
poiarization, and by considering a variety of different 
metais (M = Sc, Ti, V, Cr, Mn, Ee, Co, Ni, Cu, Zn 
and Ag). The effect of iocai eiectronic interactions on 
the centrai metai atom are aiso studied systematicaiiy. In 
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Sec. II we discuss the computational method, and present 
in Sec. Ill the magnetic properties. Then (Sec. IV) the 
transmission properties and finally (Sec. V) the spin fil¬ 
ter efficiency and electronic conductance are discussed. 
A summary is given in the concluding Sec. VI. 


II. COMPUTATIONAL METHOD 


The transport properties are investigated using non- 
equilibrium Green’s function and density functional the¬ 
ory as implemented in the SMEAGOL [HI and 
SIESTA packages [3l| . The wave functions are expanded 
in atomic orbitals with an energy cutoff of 300 Ry jH| . 
We use a double zeta plus polarization basis. Test cal¬ 
culations show no significant change when a single zeta 
basis is used for the M atom; however, when a single 
zeta basis set is used for the whole structure, the re¬ 
sults are clearly of inferior quality. The nuclei and core 
electrons are represented by Troullier-Martins pseudopo¬ 
tentials [HI- We employ the local density approxima¬ 
tion with Goulomb interaction (LDA+V). The on-site 
Goulomb interaction, I/, is applied to the d orbitals on 
the central metal ions. Eorces below 0.04 eV/A are 
achieved in the structural optimization of MPc, using 
the conjugate gradient method. Since for metallic chains 
a Peierls gap is expected around the Eermi energy, we do 
not relax the structure after attaching the leads to the 
molecule Q. We use the first two of the five-Au-atom 
leads, respectively, to calculate the surface Green’s func¬ 
tion, while the other three, closer to the molecule, are 
assumed to be part of what is usually called “scatter- 
ing region”. The Au-Au distance is chosen to be 2.89 A 
[HI, 05 [3^ 5 and the interface distance between MPc and 
Au is set to 1.37 A [H| . 

The structure of the devices is shown in Eig. [TJ As a 
general remark, we note that the Mulliken analysis shows 
a considerable charge redistribution due to the leads (Au 
atoms), mainly involving the M atom, and the N and 
G atoms close to the central atom and at the contact 
between leads and molecule. 



FIG. 1. (Color online) Structure of AuS-MPc-SAu. H, C, N, 
S, M and Au atoms are shown in cyan, yellow, grey, green, 
pink and golden, respectively. 


III. MAGNETIC MOMENTS 


The electronic and magnetic properties of isolated MPc 
and EieMPc molecules have been studied before jH| . 
where it was found that the magnetic moment (MM) 
is carried mainly by the metal atom. Below we will 
compare, in particular, isolated molecules with molecules 
connected to leads, with emphasis on the differences. In 
addition, since the central atom is a transition metal, we 
will encounter all possible d orbital occupations, which 
often give rise to considerable correlation effects. Thus 
it will be important to also investigate, not only for a 
few selected cases but systematically, the effect of a local 
Goulomb interaction on the electronic structure as well 
as on transport properties. In this context, we note that 
U often is considered to be a parameter, to be deter¬ 
mined by comparing specific physical quantities, such as 
the band gap or the magnetic moment, with experiment. 
However, the value obtained for U in such a way depends 
on the quantity considered. On the other hand, U can in 
principle be computed by constrained density functional 
theory. Both approaches appear not to be useful for the 
present systematic study where the system properties are 
strongly modified by the coupling to the leads. Thus our 
concept is to vary U systematically, within a reasonable 
range (from U = 0 to 8 eV, i.e., from weak to strong cor¬ 
relation) in order to elucidate the interaction-dependent 
trends in the changes of the system properties. Eor this 
goal, it appears sufficient to rely on the three values 0, 4, 
and 8 eV. 

The calculated magnetic moments are summarized in 
Table H We first discuss the total MMs for the AuS- 
MPc-SAu system. Eor M = Sc the MM differs slightly 
from that of the isolated molecule (1/is), and there is 
no effect of the U parameter, whereas the MM of the 
M = Ti system is higher than for the isolated molecule 
(2/is), inceasing with U. Apparently the d-d electron 
interaction of Ti enhances the charge transfer from the 
gold chain to the molecule when correlation increase. Eor 
V and Gr the MM for 1/ = 0 is smaller than that of 
the isolated molecule (3/is for VPc, 2/is for GrPc). As 
U increases from 4 to 8 eV, the MMs get closer to the 
isolated molecule value. Eor Mn junctions the MM is 
mainly located on Mn (4.8 /is for plain MnPc) and hardly 
depends on U. 

Eor the Ee case, the magnetic moment of isolated EePc 
happens to be close to 2 /is, but is found to be strongly 
increased in the junction due to the overlap between 
Au orbitals and molecular orbitals, the MM on Ee be¬ 
ing roughly doubled. Apparently this coupling induces 
a transition from a low-spin to a high-spin state, unlike 
the case of, e.g., Mn which is already in a high-spin state 
for isolated MnPc; cf. [HI, tables I and II. Thus Ee is an 
exceptional example where the coupling to the leads have 
a pronounced effect on the magnetic properties. Unfor¬ 
tunately, it seems not to be possible to explain this sur¬ 
prizing fact in terms of the standard picture of atomic d 
orbitals in a square planar crystal field (see also below). 
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TABLE L Total and metal magnetic moments (in / ib ) of AuS-MPc-SAu and AuS-FieMPc-SAu at different U 





Sc 

Ti 

V 

Cr 

Mn 

Fe 

Co 

Ni 

Cu 

Zn 

Ag 


U 

= 0 

total 

0.8 

2.9 

2.1 

2.3 

5.3 

5.9 

1.6 

1.5 

2.2 

0.4 

1.6 




metal 

0.1 

0.8 

2.9 

4.0 

4.6 

4.1 

1.3 

0.1 

0.5 

0.0 

0.2 

MPc 

U 

= 4 

total 

0.8 

3.9 

3.2 

3.7 

5.3 

5.9 

4.6 

3.2 

1.8 

0.4 

1.9 




metal 

0.0 

1.1 

3.0 

4.2 

4.8 

4.2 

2.7 

1.6 

0.5 

0.0 

0.2 


u 

= 8 

total 

0.8 

3.8 

3.2 

4.1 

5.3 

5.8 

4.4 

2.9 

1.6 

0.4 

2.1 




metal 

0.0 

1.1 

3.0 

4.2 

4.9 

4.6 

2.8 

1.7 

0.5 

0.0 

0.2 


u 

= 0 

total 

0.8 

1.6 

1.2 

4.0 

4.7 

4.1 

1.0 

0.0 

1.1 

0.0 

0.7 




metal 

0.0 

1.5 

2.4 

4.0 

4.7 

4.1 

1.1 

0.0 

0.5 

0.0 

0.2 

F isMPc 

u 

= 4 

total 

0.8 

1.5 

3.2 

4.0 

4.9 

3.4 

2.9 

1.9 

1.1 

0.0 

0.8 




metal 

0.0 

1.5 

2.6 

4.1 

5.0 

4.2 

2.6 

1.5 

0.5 

0.0 

0.3 


u 

= 8 

total 

0.8 

1.4 

3.3 

4.0 

5.0 

3.4 

3.0 

1.9 

1.1 

0.0 

0.9 




metal 

0.0 

1.3 

2.3 

4.1 

5.0 

4.3 

2.8 

1.6 

0.5 

0.0 

0.3 


For the last five columns in the table, Co ... Ag, we 
note, first of all, that the total MM of the respective 
isolated molecules are close to 1, 0, 1, 0, l/i^, respec¬ 
tively [35|. Thus it is apparent that the coupling to the 
leads has a significant effect in all cases for MPc, while 
the disturbance due to the leads is rather small for the 
fiuorinated counterparts (considering first U = 0). For 
both, MPc and fiuorinated MPc, we find a considerable 
correlation dependence for Co and Ni. However, due to 
the strong coupling to the leads for MPc, this dependence 
cannot straightforwardly explained by the electronic con¬ 
figurations of the respective isolated metal atoms. On the 
other hand, there is hardly any [/-dependence for Cu, Zn 
and Ag (filled d shells). 

Tabled] also shows that the magnetic moment on the 
central metal atom in most cases differs considerably 
from the total magnetic moment. The [/-dependence of 
the metal magnetic moment generally is small, except for 
Co and Ni. 

Because of the high electronegativity of F, the MMs of 
the AuS-FieMPc-SAu systems generally are less than or 
equal to those of AuS-MPc-SAu. For M = Cr, Cu, and 
Zn the MMs of the AuS-(Fi6MPc)-SAu junctions agree 
with those of the isolated molecules, independent of U. 
In most cases, the MMs do not change when increasing U 
from 4 to 8 eV, especially in the fiuorinated systems. The 
magnetic moments generally are slightly higher for the 
interacting case ([/ = 4 or 8 eV) than for zero Coulomb 
interaction, which is reasonable since increasing U will 
lead to a decrease of hybridization among the d electrons. 

In order to obtain some further insight into the above 
results, we also studied in some detail the charge and 
spin density isosurfaces of the considered systems. For 
example, the charge isosurfaces of isolated ScPc was pre¬ 
sented already in (Fig. 2), with the result that there 
is hardly any charge on the metal. In the transport sit¬ 
uation, this does not change; however, due to the good 
coupling to the Au leads, which is mediated by the above- 
mentioned addition of an S atom, some charge is being 
transferred to the contact region. The good coupling is 


also apparent in the spin isosurface. As a function of the 
interaction parameter [/, it is found that the spin den¬ 
sity is elongated along the transport direction, and that 
some spin density is shifted towards the leads, which is a 
general trend obtained from the SMEAGOL data. 

As another example, the spin density isosurface of 
AuS-AgPc-SAu is shown in Fig. [2j For clarity of pre¬ 
sentation, we show here and in the next figure only the 
spin isosurface for the molecule plus the first Au atom 
on both sides, i.e. we cut out the remaining part of the 
leads. Of course, the isosurfaces were computed for the 
complete system. In Fig. [2] the above-mentioned good 
coupling is clearly visible as a relatively strong spin con¬ 
tributions at and near the contact. Also apparent is the 
small magnetic moment on the metal; note the tiny con¬ 
tribution from the minority spin on Ag. Similar to ScPc, 
the U -dependence is smooth, following the general trends 
mentioned above. 



FIG. 2. Spin isosurface for AuS-AgPc-SAu as obtained from 
SIESTA (i.e., for [7 = 0), taking into account the states be¬ 
tween —2 and 0 eV; cf. Fig. [H Red: majority spin, blue: 
minority spin. Shown is only the part containing the central 
molecule and the first (left and right) Au atom. 

In Fig. [3] we present the spin density isosurfaces for 
AuS-FePc-SAu and AuS-FieFePc-SAu. Most notable is 
the fact that the spin on Fe appears to be the same for 
the two cases, in agreement with the results for the mag¬ 
netic moment, cf. Table [H For FePc, one also observes 
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large contributions to the majority spin density from N 
atoms, and from C atoms which are not bonded with 
N. However, the C atoms which are bonded with N as 
well as the central metal atom contribute to the minority 
spin density. There is hardly any contribution from the H 
atoms. Concerning the fluorinated FePc system (right), 
the spin density is clearly smaller than for FePc; cf. Table 
[H This reduction of the spin density can be attributed to 
the reduction in the spin density of C and N atoms. 



the spatial overlap region becomes larger. The changes 
in T{E) when increasing the interaction parameter are 
very small. 




FIG. 4. Left: Transmission coefficient for AuS-ScPc-SAu; 
right: corresponding projected density of states. Note that 
the An PDOS has been scaled by a factor 0.1 for easier com¬ 
parison. 


FIG. 3. Spin isosurface for AuS-FePc-SAu (left) and AuS- 
FieFePc-SAu (right; F replacing H at the lobes shown in grey) 
as obtained from SIESTA (i.e., for F = 0), taking into account 
the states between —2 and 0 eV; cf. Fig. [8l Red: majority 
spin, blue: minority spin. Shown is only the part containing 
the central molecule and the first (left and right) An atom. 


IV. TRANSMISSION PROPERTIES 

In the following, we present the transmission coeffi¬ 
cients for selected systems. Note that the d orbitals at 
the M atom have a 04 /^ symmetry with a crystal field that 
splits the orbitals as ai^ ^ 2 g i^xy)^ {^xz^ ^yz) 

and bi^ {dx 2 _y 2 ) states. For some systems, like FieScPc 
and TiPc, T{E) is found to be negligible: the projected 
density of states (PDOS) of these junctions are mostly 
formed by sharp lines, which indicates a negligible cou¬ 
pling of molecular orbitals with the atomic orbitals in the 
leads. In the following, we concentrate on a few repre¬ 
sentative cases. 

ScPc. The transmission coefficient obtained for the 
AuS-ScPc-SAu system is shown in Fig. HI The spin com¬ 
ponents are distinguished by positive and negative val¬ 
ues. A tiny shift of the T{E) peaks is noticeable as a 
function of U, see left hand side of Fig. ID The PDOS for 
U = 8 eV on the right hand side of Fig. 0] shows that the 
C and An states dominate in the energy range —2 and 2 
eV, whereas N and Sc states appear only below Ep. The 
8iig states are part of the HOMO and the eg states are 
part of HOMO —1 spin down and HOMO—2 spin up. The 
transmission peaks at —0.2, —0.3 and —0.4 eV are mainly 
due to HOMO, HOMO —1 and HOMO—2, respectively, 
which are localized on Sc and Pc. The transmission at 
— 1.3 eV is dominated by the molecular orbitals for both 
spins. The overlap between molecular orbitals and Sc 
orbitals is a signature that can be observed in T{E) as 


Fluorinated VPc. Figure [5] shows T{E) and the PDOS 
for AuS-FieVPc-SAu where all transmission peaks are 
below the Fermi energy. We find small variations in 
T{E) for energies below —0.7 eV as the value of U is 
changed. This is likely due to the fact that when U in¬ 
creases, electronic localization is favored and the trans¬ 
mission becomes more localized in a smaller energy win¬ 
dow, as shown in the left hand side of Fig. [5l Near the 
Fermi energy the direction of spin and the corresponding 
DOS strongly depend on U. The projected density of 
states for = 4 eV is depicted on the right hand side of 
Fig.O The degenerate eg states contribute to all molecu¬ 
lar states along the junction and also to the transmission 
coefficient below the Fermi energy. At —0.8 and —0.9 eV 
we find contributions mainly from the bi^ orbital. 



-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 

E-EF(eV) 



-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 


E-Ep (eV) 


FIG. 5. Left: transmission coefficient for AuS-FieVPc-SAu. 
Right: corresponding projected density of states 

We note in passing that the transmission behavior for 
AuS-MnPc-SAu and AuS-FieMnPc-SAu is similar to the 
case of ScPc and and fluorinated VPc. The transmission 
coefficient is practically independent of U, and very quite 
small. For a different lead material (namely a SWNT), 
however, it has been suggested [26| that MnPc can be 
used as spin filter. 

FePc. We now consider Fe in the center of the Pc 
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molecule. The effect of the U parameter is very clear 
in this case. Some transmission coefficient peaks happen 
to be independent of [/, in particular, those at —1.3 {eg 
and molecular states) for both spins, —0.4 eV (ai^b 2 ^ 
and molecular states), and —0.1 eV (bi^ and all other 
states of An and the molecule) for spin up. From —1 to 
—0.7 eV, some T{E) peaks appear only for 1/ = 4 eV, 
and at —0.4 eV for = 8 eV only in case of spin down. 
Furthermore, at 0.9 eV (e^bi^) but only for spin down, 
there is a single T{E) peak for = 0. For U = 0 and 
8 eV, T{E) has peaks at —0.6 eV (ei^bi^ and molecular 
states). We find zero transmission at the Fermi energy, in 
contrast to the case where FePc is connected to a SWNT 

mil. 




-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 -2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 

E-Ep(eV) E-Ep(eV) 


FIG. 6. Transmission coefficient for (left) AuS-FePc-SAu and 
AuS-CoPc-SAu (right) 

CoPc. For the case of Co in the center of the Pc 
molecule, the T{E) peaks are sharp lines, and their po¬ 
sitions vary with U. There are two peaks below and one 
peak above Ep for every U value, see right hand side of 
Fig.[ 6 l For [7 = 0, the peak in T{E) at —0.5 eV is related 
to molecular as well as An and bi^ states. There are two 
additional peaks at 0.5 eV (molecular orbitals with bi^) 
and at 1.4 eV (molecular orbitals with e^bi^ai^). 

CuPc. For the AuS-CuPc-SAu junction, several trans¬ 
mission peaks are found below the Fermi energy, and, in 
particular, they are enhanced for 1/ = 4 eV, see Fig. [3 
Above the Fermi energy, there is a single peak, which 
changes position and intensity as function of 7/, the po¬ 
sition shifting away from Ep as U increases (because of 
inceasing repulsion between electrons). The energy dif¬ 
ference between the first peak below and the first one 
above Ep for U = 4 is smaller than the corresponding 
one for U = 0 and 8 eV. 

AgPc. Last but not least we discuss AuS-AgPc-SAu 
and AuS-FieAgPc-SAu (Fig.[ 8 j). The T{E) peaks demon¬ 
strate a good coupling between lead and molecule as 
shown in the left hand side of Fig. [ 8 ] (AuS-AgPc-SAu). 
We find peaks for spin up states at Ep, independent of U, 
which indicates that this configuration has a high poten¬ 
tial as a spintronic device. The change of the interaction 
parameter does not affect the positions of the peaks be¬ 
low the Fermi energy. However, we observe a strong de¬ 
pendence of the peak position and intensity on U above 
the Fermi energy: as U increases, the peak moves further 



-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 

E-Ep (eV) 


FIG. 7. Transmission coefficient for AuS-GuPc-SAu 




E-EF(eV) E-EF(eV) 

FIG. 8. Transmission coefficient for AuS-AgPc-SAu (left) and 
AuS-FieAgPc-SAu (right) 


away from Ep, and the intensity increases almost by a 
factor of two. This indicates that as electron localization 
increases, the transmission also increases for that par¬ 
ticular state. For the AuS-FieAgPc-SAu junction, right 
hand side of Fig. [51 the transmission peak near —0.1 eV 
is finite for both spin polarizations, and independent of 
U. The transmission near the Fermi energy is dominated 
by molecular states with ai^bi^ orbitals. 

For completeness we discuss the occupation of the d 
orbitals on the metals which were addressed in detail in 
this section. First, we find only minor differences be¬ 
tween MPc and their fluorinated counterparts. Second, 
the occupation numbers of the five orbitals add up to 
numbers close to the values to be expected from the pe¬ 
riodic table: Sc: « 0.7, decreasing by about 20 % with 
increasing U; V: ~ 3.1, almost independent of U; Fe: 
« 5.5, decreasing by about 10 % with increasing U; Co: 
close to 7, independent of U; Cu: « 9.3, almost inde¬ 
pendent of 17; Ag: « 9.4, increasing by about 1 % with 
increasing U. Considering the simple—from an atomic 
point of view—case of Sc, it is already apparent that an 
interpretation in terms of the standard level picture is 
incomplete. This conclusion is further supported by the 
case of Fe, which has a rather high magnetic moment as 
discussed above: we find that all five d orbitals are oc¬ 
cupied by close to one up electron, while the spin-down 
contribution is mainly due to an Og orbital (« 0.5 for 
77 = 4 eV, dropping to ~ 0.1 when U is increased to 8 
eV, consistent with the increase of the metal magnetic 
moment from 4.2 to 4.6cf. Table [I]). For the V sys- 
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tern, the d orbital occupations depend only weakly on 
interaction: the eg (up) and bi^ (up) occupations are 
slightly increased when U is increased, while the others 
decrease by a small amount. For the Co system, we find 
the up occupations to be close to one, while the occu¬ 
pations of the h 2 g (down), ai^ (down), and one of the 
eg (down) orbitals are considerably smaller. For Cu and 
Ag, all occupations are close to one, except for the ai^ 
(down) orbital (^ 0.6 and ^ 0.7, respectively). 


V. SPIN FILTER EFFICIENCY AND 
CONDUCTANCE 


The spin filter efficiency (SFE) for an electronic de¬ 
vice, defined as 


mEF)+T^{EFy 


( 1 ) 


is the ability of a device to pass a particular spin compo¬ 
nent. Here, T^{Ef) and Ti{Ef) denote the transmission 
coefficients of spin up and spin down electrons at the 
Fermi energy, respectively. Our findings are summarized 
in Table im the positive values corresponding to the case 
where “up” is the majority spin, while a negative sign 
means that the majority spin is “down”. We note that 
the magnitude and the sign of the spin filter efficiency 
can vary with the interaction parameter. In some cases 
the majority spin in FieMPc is opposite to the majority 
spin in MPc, e.g., for Sc and Cu. Our results suggest 
that junctions containing Sc and Ag, as well as FePc, 
ZnPc, FieTiPc, FieVPc and FieCuPc could potentially 
be used in spintronic devices. Note, however, that ScPc is 
chemically unstable (sH, • In several of the considered 
molecules we observe zero spin filter efficiency (indepen¬ 
dent of U), in particular, for CrPc, CoPc, NiPc, and their 
fluormated counterparts. 

However, a reasonable electronic conductance is also 
essential for the performance of a device, hence we deter¬ 
mine in addition the conductance, 

G = Go iFiEF) + FiEF)], Go = 2eyh. ( 2 ) 


For the uncorrelated case, 1/ = 0, the junction which in¬ 
cludes AgPc shows the highest conductance, close to 0.22 
Go- The next higher one is CuPc with 1.8X10“^ Go- This 
value is smaller than the experimental result found in aro¬ 
matic molecules by two orders of magnitude [ssl, [39| , but 
our result roughly agrees with another theoretical cal¬ 
culation [25|. Differences in theoretical results generally 
depends on the construction of the junctions, as well as 
the employed calculational method and the details of the 
(weak) interaction between molecule and leads. From a 
practical point of view, the results will also depend on 
the interaction between the junction and the substrate, 
which is not taken into account here and in most other 
studies. 

The conductance of ZnPc is found to be 3 x 10“^ Go, 
which is the same order as previous experimental and 


theoretical results for an oligo-porphyrin molecular wire 
by using Zn in the molecule center [^. We mention in 
passing that for TiPc and VPc, which are also chemically 
unstable [3^ . [4Q| - [4^ . a vanishing conductance is obtained. 

In FieMPc, the conductances in most cases are less 
than for the corresponding MPc. The highest conduc¬ 
tance is found for FioAgPc, namely 0.002 Go, the next 
highest value being 1.5 x 10“^ Go for FioFePc. In addi¬ 
tion, the conductance of FioNiPc and FieCuPc is of the 
same order as the conductance of FioFePc. 

The electronic conductance is somewhat insensitive to 
an increase of U for all systems, in particular, it remains 
zero for thoses cases where we find zero conductance at 
U = 0. On the other hand, when the conductance is 
finite at 1/ = 0, we find it either to remain constant or to 
decrease with increasing U. For example, for 1/ = 4, the 
conductances for CuPc, ZnPc and AgPc are 6 x 10“^, 3 x 
10“^ and 1.1 X 10 ^Go, respectively. When we increase 
U to 8 eV the conductances for the same systems drop 
to 0, 3 X 10“^ and 7.3 x 10“^Go, respectively. 


VI. SUMMARY 

In this paper we studied systematically the magnetic, 
electronic and transport properties of metal phtahlocya- 
nines and fluorinated metal phthalocyanines connected to 
An leads, including some MPc molecules which are chem¬ 
ically unstable. We employed the LDA+U approach, 
with values of the Coulomb interaction parameter rang¬ 
ing from U = 0 to U = 8 eV. The magnetic moments 
are largely determined by the hybridization between d 
metal and An states near the Fermi energy. The mag¬ 
netic moments (Table [J) to some extent vary with the 
electron-electron interaction on the central metal atom, 
U, which also for some systems considerably modifies the 
spin filter efficiency (Table [II]). Considering Eq. ([T]), it is 
apparent that the spin filter efficiency is a very sensitive 
quantity whenever the transmission coefficients for both 
channels are very small. 

In particular, the magnetic moments are found to in¬ 
crease with increasing correlation, or are roughly inde¬ 
pendent of U. For the transport properties, the situation 
is less clear: it appears that for some systems the spin fil¬ 
ter efficiency hardly depends on U, while others are quite 
sensitive to correlation effects. There can be a remarkable 
difference between MPc and its ffuormated counterpart. 
In detail, the results can be explained by the respective 
electronic structure near the Fermi surface, which, how¬ 
ever, is not only determined by the orbitals of the central 
metal atom but also by the coupling to the leads. In order 
to finally clarify the role of electronic correlation in the 
series of junctions considered, additional experimental ef¬ 
forts will be needed. Overall our theoretical results com¬ 
pare well with other theoretical and experimental results 
for the electronic conductance where these are available. 

The structures with Ag in the center of the junction 
are a notable exception: they show both a good spin filter 
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TABLE IL Spin filter efficiency for AuS-MPc-SAu and AuS-FieMPc-SAu in % 




Sc 

Ti 

V 

Cr 

Mn 

Fe 

Co 

Ni 

Cu 

Zn 

Ag 


U=0 

-18.0 

0 

0 

0 

0 

99.9 

0 

0 

0 

0 

99.9 

AuMPc 

U=4 

-3.0 

0 

0 

0 

0 

51.7 

0 

0 

8.3 

68.4 

99.9 


U=S 

-16.0 

0 

0 

0 

0 

-2.5 

0 

0 

0 

72.9 

97.9 


u=o 

48.8 

99.9 

14.2 

0 

0 

0 

0 

0 

0 

0 

99.5 

AuFieMPc 

U=4 

70.4 

98.0 

64.3 

0 

0 

0 

0 

0 

-99.9 

0 

99.2 


u=s 

74.1 

99.0 

-90.6 

0 

0 

0 

0 

0 

-99.3 

0 

98.6 


efficiency and a reasonable electronic conductance since 
the contribution of Ag states at the Fermi energy is large. 
Hence AgPc and FigAgPc junctions can potentially be 
used in spintronic devices. 
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